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Abstract 


The effects of pion (n) production are expected to play an impor- 
tant role in radiation exposures in the upper atmosphere or on the 
Martian surface. Nuclear databases for describing pion production 
are developed for radiation transport codes to support these studies. 
We analyze the secondary energy spectrum of pions produced in 
nucleon-nucleon (NN) collisions in the relativistic one-pion exchange 
model. Parametric formulas of the isospin cross sections for one-pion 
production channels are discussed and are used to renormalize the 
model spectrum. Energy spectra for the deuteron related channels 
(NN — » dn) are also described. 


1. Introduction 

A convenient representation of the differential cross section in energy of particles created in the 
interactions of space and atmospheric radiation with materials is required for radiation transport corn- 

reactions, an abundance of mesons is produced through the nuclear force. High-energy reactions 
increase in relative importance in the upper atmosphere due to the Earth’s magnetic field, reducing 
lower energy ion components which have insufficient energy to overcome the threshold for meson pro- 
duction. The one-pion production channels dominate the meson production for energies up to about 
1.5GeV/amu and represent an important contribution at higher energies. The energy region below 
1.5 GeV/amu extends above the galactic cosmic ray (GCR) peak at 0.2 to 0.6 GeV/amu. Also, most of 
the energy range of particles seen in solar particle events (SPE) are dominated by the one-pion exchange 
interactions. The scattering of high-energy protons or neutrons on hydrogen is important because pro- 
tons and neutrons represent over 90 percent of the particle flux in materials and also because hydroge- 
nous materials in tissue and material structures are important. The description of the one-pion 
production mechanism in nucleon-nucleon (NN) collisions is also needed for modeling the pion produc- 
tion mechanism in inclusive proton and neutron collisions on target atoms in applying reaction theory or 
in Monte Carlo simulations (refs. 4 and 5). 

In this report, we describe the calculation of the secondary energy spectrum in proton-proton (pp) 
and neutron-proton (np) reactions. The one-pion production channels are modeled by using the relativis- 
tic one-pion exchange (OPE) model. Parametric models for the isospin components in NN reactions and 
jcN reactions are discussed. The one-pion channels discussed here can be appended with parameteriza- 
tions of high-energy models (refs. 6 through 8) to provide the energy spectrum for inclusive pion pro- 
duction in NN collisions at overall energies of interest for space radiation studies. 


2. One-Pion Production Cross Sections 


The cross sections for production of a single pion in nucleon-nucleon (NN) collisions may be 
written in terms of four independent cross sections by applying isospin conservation. For a transition 
from an initial isospin state / to a final isospin state Ij the cross section is denoted a/./ . For formation 
of a deuteron d in the final state, a superscript d is used. In table 1 we list the isospin and masses of the 

the various reaction channels in NN collisions. Various authors have considered parameterizations of 
these cross sections: a 10 , a 01 , a n , and cf (J . The fits of Wilson and Chun (1988) are useful because 
they extend over all energies. The work of VerWest and Arndt (ref. 9) is significant because they have 




reconsidered discrepancies in older data sets, resulting in a large reduction in the isospin zero cross 
section a 01 over previous estimates and are in good agreement with the more recent measurements 
(ref. 9). The formula of VerWest and Arndt (ref. 9) for o^ 0 is (in units where 7i = c- 1) 
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with the nucleon mass, m N = 0.939 GeV, m n the pion mass, and m K = 0.138 GeV. The formula used 

in reference 9 for a , , is 
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where 


Z + = (,js-m N -M o )(2ir 0 ) 

Z_ = (m N + m K -M 0 )(2/r o ) 

with M g = 1.22 GeV and V o = 0.12 GeV. The parameters of the model are listed in table 3. The for- 
mulas of equations (2.1) and (2.2) are valid only for T L < 1.5 GeV. 

Wilson and Chun have considered the following parameterizations of the isospin cross sections for 
all energies: 
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where T th is the single-pion production threshold and T L is the laboratory nucleon energy. 

Comparisons of the fits described above are shown in figure 1. Large discrepancies exist, especially 
for o Q1 and c 10 . Because the model of VerWest and Arndt (ref. 9) is more accurate at lower energies, 
we will use this model below 1.3 GeV. The resulting fit is shown in figure 2, and comparisons to data 
(ref. 10) for K + and n° production in proton-proton (pp) collisions are shown in figure 3. The result of 
joining equations (2.1) and (2.2) with equations (2.3) to (2.6) will represent the one-pion production 
data quite accurately over the energy range of interest for radiation transport codes and will be used to 
renormalize the model spectra described in the next section. 


3. Pion Energy Spectrum 

The deuteron production channels are two-body final states and are therefore much easier to param- 
eterize than the other pion production channels. By isospin conservation (neglecting Coulomb effects), 
we have the relationship 
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The angular distribution is parameterized as 
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where N^q is a normalization constant and 0 is the cm scattering angle. The energy dependent parame- 
ters A% and are 
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and 
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Comparisons of equations (3.1) to (3.4) to experimental data (refs. 11 and 12) are shown in figure 4(a) 
for the K + d — > pp reaction and in figure 4(b) for the pp —> n d reaction. 

The energy distribution in the laboratory frame of the final deuteron is related to the center of mass 
(c.m.) angular distribution by 


do 

dT d 



(3.5) 


where the function X is defined 


X(s,A,B) = ( s-A-B ) -AAB 

2 


(3.6) 


and s is the Mandelstam variable given by s = + P/v 2 ) ■ 

In applying equation (3.5), we use 
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where barred quantities are the c.m. values given by 
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The kinematical limits on the kinetic energy of the deuteron T d are found from equation (3.7) by observ- 
ing that | cos 9 < 1. 


The energy distribution in the laboratory reference frame for the pion is 
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and a similar expression is found for the deuteron spectrum. Calculations of energy spectra of secondary 
pions and deuterons for several beam energies are shown in figures 5 and 6, respectively. 

The NN — » NNtt reactions are assumed to proceed through the formation and decay of the A reso- 
nance. The mechanism is described by figure 7(a). The A forms an isospin quartet (7 = 3/2), and it is use- 
ful to consider the coupling of the kN system wave functions in isospin space to understand the A decay 
properties. Denoting the nucleon isospin by In and the pions by I n , we have 

IMO = £l"z><"zl W-V < 3 - 13 > 

A/ z 

Equation (3.13) is used to obtain the components of the A wave function with the result 
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(3.14c) 
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Branching ratios for the formation of various pion species in NN collisions 
equations (3.14). 

are obtained by using 

The invariant differential cross section for NN — » NA has been evaluated in the relativistic one- 
pion exchange (OPE) model. We use this model and the assumption of isotropic decay of the A in its 
rest frame to obtain the momentum distribution and energy spectrum of pions and nucleons in the 
NN — » NNtc channels. The Mandelstam variables for the reaction are defined in terms of the four 
momentum vectors of the various particles defined in figure 7 and given by 

^ = (K l + K 2 ) 2 

(3.15) 
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The cross-section distribution in t is written in terms of the matrix element M as 


do 

dt 
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The matrix element M is decomposed as 

!WI 2 = |M dir | 2 + iM CI | 2 + lM INTi 2 


The direct amplitude is given by (ref. 13) 



(3.20) 


(3.21) 


where / and /* are the coupling constants with values 1 .008 and 2.202, respectively, and m n and 
m A are the pion-fixed A mass with values 0.139 GeV and 1.232 GeV, respectively. In equation (3.21), 
F(t) is the form factor for the off-shell meson which is parameterized as 


6 




F{t) = 


(3.22) 



where the value of the parameter A will be discussed below. The exchange term in equation (3.20) is 
equivalent in form to equation (3.21) with the replacement (u <-> t). 

The interference term is given by 



(3.23) 


In order to account for the finite width of the A resonance, a mass distribution p(p ) is introduced 
through 


= A(s,r, M )p(p 2 ) (3.24) 

dt d\i 


2 2 

where p = (AT 4 + K n ) and A(s, t, u) is the invariant cross section of equation (3.18) with the fixed 
A mass m A replaced by p. The mass distribution is parameterized in terms of the elastic pion-nucleon 
cross section and the A width as (ref. 1 3) 
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The width is parameterized as 
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with the form factor Z given by 


Z(V 2 ,ml) = 



(3.28) 


and parameter values V o = 0.12 GeV, and k = 0.2 GeV. By neglecting Coulomb effects, the shape 
of the secondary energy spectrum for the components of the isospin triplet of pions is largely deter- 
mined by the nN cross section in equation (3.25), with the remaining kinematic factors in 
equations (3.18) to (3.28) being identical for each pion species. The description of the kN cross 
sections is given in the appendix. 


The decay of the A (A —> Nn) is assumed to be isotropic in the A rest frame such that the pion 
distribution from the decay is 


do 

dt d\l 2 d£l K 


t,u)p(\l 2 ) 


and the nucleon distribution from the decay is 


do 

2 * 

dt d\x dtl N 


^A(s, t,«)p(p 2 ) 
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(3.30) 


where starred quantities are in the A rest frame. A second contribution to the nucleon spectrum comes 
from the nondecaying nucleon line in figure 7, as discussed below. The A mass is found as 
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allows the invariant momentum distribution of the pion to be written 


E « 7T = —zl ^ 

d Pn An 1 


P\P7> 


W2 


— A(s, t, «)p(p 2 ) 


^ Pn > 


BE, 


(3.33) 


8 




with a similar expression for the decay nucleon spectrum. The energy spectrum of the pion in the labo- 
ratory system is then 


do 
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= —2 1 da l 


n 4n 
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The laboratory energy spectrum of the decay nucleon is 
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The energy spectrum of the recoil nucleon is given by (ref. 14) 
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where p a is the laboratory momentum of the incident nucleon. 

A comparison of equation (3.24) to experimental data for the pp — » nA reaction (refs. 10 and 15) 

is shown in figure 8, with good agreement found. The parameter A in the form factor of equation (3.22) 
has a strong effect on the calculations. As noted by Jain and Santra (ref. 14), the inclusion of distorted 
waves effects, the value of A chosen, and a value of — 1 GeV provides the best fit to data. Calculations of 
energy spectrum for ti + and n° production in pp collisions and 7t + , K°, n~ production in neutron-proton 
(np) collisions are shown in figures 9 and 10, respectively. The variations in the A — > nN decay vertices 
due to the isospin dependence provide only a modest change in shape between the various production 
channels. Coulomb effects which are not treated herein will provide further dependence on the pion 
charge. 

4. Inclusive Pion Production Spectrum 

As the kinetic energy of nucleons increases, the threshold for 2, 3 . . . pion production is reached and 
so is production of heavier mesons such as the kaon. The threshold energies for several production pro- 
cesses are listed in table 4. The production threshold is also dependent on charge conservation when 
individual species of mesons in the reaction are produced. Inclusive meson production data have been 
parameterized in convenient forms by several authors (refs. 6 through 8). The spectrum in one-pion, 
two-pion, and other channels will be somewhat distinct due to the mechanisms involved. The two-pion 
production channels have been considered by Stemheimer and Lihdenbaum (ref. 16) using a purely 
kinematic form of the isobar model. There are two distinct mechanisms for two-pion production which 
are through the excitation of two A’s: 

N + A— » Aj + A 2 ->N + yV + jt + 7t (4.1) 


or through the excitation of higher mass nucleon resonances: 


N + N^>N + N*-+N + N + 7t + Jt 


(4.2) 
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where N* is a nucleon resonance of higher mass than the A, which is assumed to decay through the 
emission of two pions. The mechanisms of equations (4.1) and (4.2) are illustrated in figures 7(b) 
and 7(c), respectively. The mechanism of equation (4.2) will have a cross-section structure similar to 
the one-pion model, NN NNrc , described above. However, with A — > Nk, the vertex is replaced by 
the N* — » Nun vertex. The mechanism of equation (4.3) will contain the mass distribution of two A’s 
with a structure such as 


2 = < 4 - 3 ) 

dt d\i { d\i 2 v y V 7 

thus requiring one additional numerical integral to obtain the pion energy spectrum. 

The inclusive pion spectrum can be represented as a sum over the spectrum for each multiplicity at 
pion as 


da_ _ da_ t/a 
dE n - dE ni dE^ 


(4.4) 


where 


da 

dE„ 


represents the spectrum described by equation (3.34). The inclusive spectrum can be repre- 


sented by the one-pion channel plus the high-energy model of Schneider, Norbury, and Cucinotta 
(ref. 7) renormalized to exclude the one-pion contributions, or alternatively, by including the two-pion 
production channels separately, as described previously. 


5. Concluding Remarks 

One-pion production channels will dominate pion production for a significant fraction of galactic 
cosmic ray (GCR) and solar particle event (SPE) exposures in free space, in the upper atmosphere, or on 
the Martian surface. In this report, the one-pion production cross sections were discussed, and a conve- 
nient formula for their numerical representation was found. These models will be appended with high- 
energy models to span all energies of importance in GCR studies. The formula described here also can 
be used to model pion and nucleon production spectra in nucleon-nucleus and nucleus-nucleus 
reactions. 

NASA Langley Research Center 
Hampton, VA 23681-2199 
September 1, 1998 
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Appendix 


Pion-Nucleon Scattering 

The nN scattering amplitude is written to represent the pion isospin I n and nucleon isospin as 

f nN = §(/ 1/2 + 2/3/2) + ■ X N ( A1 ) 

where /j /2 and fy 2 are amplitudes for total isospin 1/2 or 3/2, respectively. Introducing the total isospin 
of the 71 N system, 

T = I n+ l -z N (A2) 

and nN wave functions 

I TT Z ) = £<Vz% x Z I TT Z> I Vz x /V x Z> ( A3 > 

and noting 

< A4 > 

leads to the following relations for elastic scattering: 


<* + P 1 f nN 1 P) = f 3/2 

(A5) 

(n°p | f KN | K°p) = \{f 1/2 + 2f 3/2 ) 

(A6) 

( n P 1 f nN 1 n P) = g(/ 3/2 + 1 / 2 ) 

(A7) 

(7T + « I | 7T + n) = |(/ 3/2 + 2/ 1/2 ) 

(A8) 

<71°« | / KW 1 7i°n> = |(/i/2 + 2 / 3 / 2 ) 

(A9) 

(n n | f nN | n n) = / 3/2 

(A10) 


with similar relations found for charge-exchange matrix elements. 
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The total cross sections can be found from the optical theorem. Numerous measurements exist for 
the total 7t p and n p reactions. From equations (A5) to (A 10), we can find solutions for the other 
hN collision pairs (neglecting Coulomb effects) by using equations (All) through (A 14). 


71 p 
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CT. "f 0 _ 0ry 

K p Up EX 


(tt p) 


(All) 


T T 

ct + = a _ (A 12) 

71/1 n p 

T T 
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T T 

a _ = o + (A 14) 

n n n p 
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Table 1 . Isospins and Masses of Particles 


Particle 

Isospin 

Isospin projection 

Mass, GeV 

P 

1/2 

1/2 

0.9383 

n 

1/2 

-Ml 

.9396 


1 

1 

.1396 

n° 

1 

0 

.1350 

K 

1 

-1 

.1396 

A ++ 

3/2 

3/2 

1.232 

A + 

3/2 

Ml 

1.232 

A" 

3/2 

-Ml 

1.232 

A' 

3/2 

-3/2 

1.232 

d 

0 

0 

1.876 


Table 2. Isospin Components of One-Pion Production Channels in 
Nucleon-Nucleon (NN) Collisions 


Channel 

Isospin cross section a, 7 

P p 

pp — > n + d 

d 

CT I0 

pp -» ppn° 


+ 

pp — > pnn 

a,i + a l0 

o j 

np — > 7t d 

1 d 

2 °io 

np — > npn° 

1. x 

2^ a io + a oi) 

np — » nnn 

i(CTii + CF 01 ) 

np — » ppn 

i(CTil + CT 0 l) 


Table 3. Isospin Cross-Section Parameters 
[From ref. 9] 


Parameter 


„ MeV 
MeV . 


„d 

°10 

°u 

a 10 

°01 

6.030 

3.772 

15.28 

146.3 

1.700 

1.262 

0 

0 

1203 

1188 

1245 

1472 

134.3 

99.02 

137.4 

26.49 


Table 4. Meson Production Thresholds 

Channel Threshold, GeV 

NN -> NNtc 0.29 

NN -> NNtctc .602 

NN -> NN7C717C .934 

NN — » KA 1.12 
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Figure 1 . Comparison of isotopic cross sections o for one-pion production versus laboratory kinetic energy of 
incident nucleon. 
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Figure 2. Comparison of isotopic cross sections a for one-pion production versus laboratory kinetic energy of 
incident nucleon using connection of low- and high-energy fit equations. 
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collisions versus laboratory nucleon kinetic energy. 
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T n = 188 (Tp = 662 MeV) 
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(a) 7t+ + d — » p + p reactions. 
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(b) p + p — » + d reactions. 

Figure 4. Comparisons of model fits to experimental data (ref. 1 1) for angular distributions; T n is kinetic energy of 
pion; T p is kinetic energy of proton. 
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Figure 6. Secondary deuteron d spectrum in p + p -> n + + d reactions for several proton beams. 
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Figure 7. Diagrams for one-pion production through A or nucleon resonance /V* formation and decay in nucleon- 
nucleon (NN) collisions. 
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Figure 8. Comparison of model to experiments (refs. 10 and 15) for invariant distribution in pp -> n& 


reactions. 
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Figure 10. Calculations of pion (one-pion contribution) energy spectrum in neutron-proton (np) collisions at 
several beam energies. 
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